Light-meson properties from the Bethe-Salpeter equation 
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Abstract. We discuss how to extract observables from an inhomogeneous vertex Bethe-Salpeter equation without resorting 
to the corresponding homogeneous equation. As an example we present a prediction for the e + e~ decay width of the p(1450) 
or p' meson. We also attempt to identify the momentum range contributing to a vector meson's decay constant. 
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INTRODUCTION 

Due to their quantum numbers, vector-meson resonances 
are easily produced in e + e~ scattering. In addition, the 
p is among the lightest mesons and thus object of nu- 
merous studies and a prime target for theoretical in- 
vestigation. In QCD, the Dyson-Schwinger-equation ap- 
proach offers a nonperturbative continuum method to 
study mesons as bound states of quarks and gluons via 
the Bethe-Salpeter equation. Herein we present some 
new results in the light of recent progress regarding the 
methods used to approach such a bound-state problem. 

THE VERTEX BSE 

A general vertex F(q,P), the inhomogeneous Bethe- 
Salpeter amplitude (iBSA), that connects quark and anti- 
quark to a color-singlet current satisfies the equation 

F(p,P)=r + j —^K(P,p,q)S a (q + )r(q,P)S b (q-) 

(1) 

(inhomogeneous or vertex BSE), where K is the quark- 
antiquark scattering kernel, and To a current which de- 
fines the channel under investigation by spin, parity 
and charge-conjugation parity. In this work, we investi- 
gate vector quantum numbers, J c = 1 , such that we 
choose (cf. [1]) 

r =3iy M , (2) 

where the renormalization constant Zi is calculated from 
the gap equation, cf. [2], 

The iBSA has poles at the positions of the respective 
bound states, and it can be decomposed as 

{q,P) = £ ,A + regular terms , (3) 

i i 

where Tiu is the homogeneous BSA of the meson under 
investigation, Pf = —Mf the total momentum squared of 



the excitation i in the respective channel, and .ylj denotes 
a normalization constant. 

We work in the well-established setup of the rainbow- 
ladder truncation and the effective interaction proposed 
by Maris and Tandy [3], with light-quark masses and the 
parameter co = 0.4 GeV (except where noted), as given 
in [4]. 

MASSES AND DECAY CONSTANTS 

According to Eq. (3), each bound state results in a pole 
in the iBSA F^(P,q). As described in [5], the amplitude 
is decomposed into components Fi(P 2 ,q 2 ,P ■ q) and co- 
variants T^(P,q,y) according to 

r^p 2 ,q 2 ,p- q ) = Y. F i( p2 'i 2 ' p -i) T ^^r) > ( 4 ) 

i 

where the standard vector covariants [4] are orthonor- 
malized according to Tr[7^ [P, q,y)T^{P,q,y)] = Sij. It is 
well-known how to obtain meson masses and decay con- 
stants using the corresponding homogeneous BSE. Here, 
however, we explore a different approach. 

We calculate the bound state masses by fitting the 
inverse of the first component F\{P 2 ,0,0), as shown in 
Fig. 1. We obtain 

m p = 0.774 GeV m p , = 1 .034 GeV (5) 

which is in agreement with [4] and [6]. 

The decay constant for a vector meson is given by [3] 

Mif v {Pf) = y / ^Tr^fo+rj, , 

(6) 

where the trace runs over color and Dirac-indices. To 
extract the same information from the iBSA, we first 
consider the general case of a projection f-^ (P 2 ) of the 
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FIGURE 1. The inverse of the first component of the inho- 
mogeneous vector amplitude l/F\ (P 2 ,0, 0), as a function of the 
square of the total momentum P 2 . The zero-crossings give the 
masses of the ground state and the first excitation. 



iBSA r(q,P) on a current f , 



(2*)« 



The poles in F(q,P) translate into poles of f~! h \p 2 ). In 
order to calculate the corresponding on-shell projection 
ff(P 2 ) (which ultimately gives the decay constant), the 
inhomogeneous BSE is rewritten as [1] 







| 34534 


- 






V 34532 








34530 




-1.5 - 


1.0 


-0.5 


0.5 1.0 



FIGURE 2. The projection f^ h] (P 2 ) defined in Eq. (7). The 
line represents a pole fit to the data points, where the pole 
positions were taken from the fit of Fig. 1 . 



In the case of f,{Pf) the inhomogeneous term and the 
current are identical, Tq = T = ^ijp- The decay constants 
can therefore be calculated using only the residues r, of 

fy^ (P 2 )> which are extracted from the pole fit shown in 
Fig. 2, where the masses (5) are used as input. Compar- 
ing Eqs. (10) and (6), we find 



MP?) 



3(-P 2 



(11) 



r(p,p) = r - 



■^M(P,p,q)S a (q + )T () S h {q-)- 



(8) 

M(P,p,q) denotes the fully amputated quark-antiquark 
scattering matrix which contains the bound state poles 
and may therefore be written as [1,7] 

M(P,p,q) = £ pT^k + re S- temS > 

i i 

(9) 

with the index i labeling the bound states in this channel, 
and r^j denoting the charge conjugation of r^j (see, e.g., 
[2]). Inserting Eqs. (8) and (9) in Eq. (7), we find 



I / d - 1 ^K[rSa{p + )r [h] {pA)s b (p-)] 



{lit) 

TrflA] -Pi)Sg (q+)T S b (q- ) 
p2_p2 



reg. terms . (10) 



Thus, in order to obtain ff{Pf), one has to divide the 

residue of f~! h ' (P 2 ) by the square root of the correspond- 
ing residue of the projection on the inhomogeneous term, 

/r?V)- 



PREDICTION FOR f p , 

The resulting decay constants of the p -meson and its first 
radial excitation corresponding to the masses in (5) are 



fp =0.213 GeV f p , = 0.079 GeV 



(12) 



which perfectly agrees with the results from the corre- 
sponding homogeneous BSE. To arrive at a concrete pre- 
diction for the decay width of the p' into e + e~ we need 
to investigate the sensitivity of these results to the char- 
acteristic parameter of the model of Ref. [3]. Indeed, the 
excited-state result exhibits a considerable dependence, 
while for the ground state the results for both m and / 
are rather solid, i.e., a dependence is observed but small 
(of the order of 2% for m and 7% for /, see [4, 8, 9] for 
in-depth information and a thorough discussion) 

Here the important point is that one can use the sys- 
tematic behavior of this dependence to determine both a 
preferred parameter value or range for an ideal descrip- 
tion in a phenomenological sense as well as to provide an 
estimate of the systematic effects in the calculation. As 
'best parameter' in our case we extract co = 0.5 GeV and 
obtain 



m p = 0.762 GeV /„ =0.218 GeV 



(13) 



for the ground state, which translates into a width T e + e - 
of 6.92 keV. The PDG [10] quotes an experimental value 
of 7.04 ± 0.06 keV. For the p' the experimental situation 
is less clear: The PDG quote but don't use two results of 
the order of and smaller than 0.1 keV. Our excited-state 
result corresponding to (13) lies in a range accessible to 
us only via extrapolation techniques due to the analytic 
structure of the quark propagator (see, e.g., [11] for 
a discussion). However, we use extrapolations both in 
model-parameter and momentum space to reduce the 
uncertainty. We find for p' 



dfXP 1 ) 



fp, = 0.095 ± 0.039 GeV T e + e 



0.94 ± 0.66 keV . 

(14) 



DISASSEMBLING f p 

To obtain further information on the structure of the 
p-meson and its radial excitation, the contributions to 
the decay constants from the different momentum scales 
are investigated. Therefore, we solve in addition to the 
inhomogeneous also the homogeneous BSE, from which 
the decay constant can be calculated according to Eq. (6), 
and define the density df v (p 2 ) via 



UPf)=(\/^P?)jdp 2 df r {p 2 ) 



(15) 



df r (p 2 ) is plotted in Fig. 3 for the ground state and first 
radial excitation of the p-meson. The main contribu- 
tions are centered in the mid-momentum regime around 
p 2 =0.1 GeV 2 , and neither the UV nor the IR have a 
strong influence. This supports [12], where it was shown 
that fp (among other quantities) is insensitive to the be- 
havior of the effective interaction in the far infrared. The 
same is true for the (perturbatively determined) UV do- 
main. Therefore, an effective interaction, which neglects 
this, quite naturally gives reliable results for the decay 
constant as well, cf. [13]. 

We further note that the scale dependence of anal- 
ogous projections has also been investigated in lattice 
QCD [14]. 



CONCLUSIONS 

We have demonstrated how the masses and decay con- 
stants of the p meson and its radial excitation can be cal- 
culated from the inhomogeneous (vertex) Bethe-Salpeter 
equation without using any information from the corre- 
sponding homogeneous solutions, which is in contrast to 
a previous study [15]. We have exemplified our method 
via a reasonable result for f p i . 

Furthermore, we have investigated the contributions to 
fp and f p i from different momentum scales, and found 
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FIGURE 3. Plot of the density df v (p z ) for the vector ground 
and excited states. 



significant contributions from neither the far IR nor the 
UV domain. 
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